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ABSTRACT 

The classification of galaxy mergers and isolated disks is key for understanding the relative impor¬ 
tance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic 
properties of galaxies as traced by emission lines have been used to suggest the existence of a sig¬ 
nificant population of high—z star-forming galaxies consistent with isolated rotating disks. However, 
recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. 
To further investigate the robustness of merger/disk classifications based on kinematic properties, we 
carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: 
from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially 
redshift the Wide Field Spectrosgraph (WiFeS) observations of these local (U)LIRGs to z = 1.5 to 
make a realistic comparison with observations at high—z, and also to ensure that all galaxies have the 
same spatial sampling of ~ 900 pc. Using both kinemetry-based classifications and visual inspection, 
we find that the reliability of kinematic classification shows a strong trend with the interaction stage 
of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties com¬ 
pared to isolated disks, whereas a significant population of the interacting disks and merger remnants 
are indistinguishable from isolated disks based on the measurements of kinematic asymmetries or 
visual inspection. The high fraction of late-stage mergers displaying disk-like kinematics reflects the 
complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified 
disks and mergers depend on the definition of kinematic asymmetries and the merger/disk classifica¬ 
tion threshold when using kinemetry-based classifications. Our results suggest that additional merger 
indicators such as morphological properties traced by stars or molecular gas are required to further 
constrain the merger/disk classifications at high—z. 

Subject headings: galaxies: kinematics and dynamics—galaxies: structure—galaxies: interactions 


1. INTRODUCTION 

Galaxy interactions play a key role in driving the 


(e.g. 

Toomre & Toomre 1972 Toomre 19771 White & 

Rees 

19781 ISchweizer |1982| |Harnes||1988| |Wright et al. 

1990 

i3arnes|1992p. During th 

3 merger processes, gas in 
momentum and falls to- 
nse starbursts and active 

disk galaxies loses its angular 
ward the center, inducing inte 
galactic nuclei (AGN) activities 

(e.g. Barnes & Hernquist 

1996 

Mihos & Hernquist 1996 

banders &: Mirabei|1996p. 


then disperse the surrounding gas, halt t he black hole 


growth, and quench the sta r formation (e.g. Sanders et al. 
1988[ Hopkins et al.||2006 ). 

One essential way fo identify interacting galax- 


optical morphological properties (e.g. 

Surace et al.||1998 

Farrah et al.|2001 

Veilleux et al.|2002 

. IN umerical simu- 


ral galaxies can result in various observational fo otprints 


of interactions such as bridges and tidal tails (e.g. Toomre 


& Toomre|1972 Barnes & Hernquist|1992 ). While visual 
classihcation is ^ ill used tor identifying mergers out to 
1 — 3 (e.g. Be]let_alj2005 Dawra_e^^y20^ Kartal- 


tepe et al. 2010 ' Hung et al.|2U13 ' Kartaltepe et al.|2U12 ) 
various automatic classihcation schemes have also been 
developed and applied to large extraga l actic surveys (e.g. 
Abraham et al.| |2003[ IGonselice |2003[ |Lotz et ari|2004[ 


2013p . Identifying such 


Law et al. 2007 Freeman et al. 
merger signatures with either visual or automatic classi¬ 
fications at z ~ 1 — 3 is challenging du e to surface bright¬ 


ness dimming and band -shifting (e.g. Hibbard & Vacca 
1997 Petty et al. |2009 ), which can lead to a sigmhcant 


underestimation of the occurrence of mergers (e.g. Abra¬ 
ham et al. 119961 [Overzier et al.||2010[ [Hung et "S 2014[ 
Petty et al. 2014^7 Meanwhile, z ~ 1 — 3 star-forming 


galaxies may show more irregular and clumpy structure 
compared to local spiral galaxies (Elmegreen et al.[2004 


2007 

sihcation of mergers an^ 


Dekel et al. 2009^, further complicating the clas- 
d isolated disks. 


On the other hand, the kinematic properties of galaxies 
traced by ionized gas or molecular gas have the poten¬ 
tial to unambiguously distinguish mergers from isolated 
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disks (e.g. Swinbank et al. 2006 Tacconi et al. 2006_ 
Menendez-Delmestre et al.| 2013|). integral tieid spec- 
trograph (ii'S) surveys witn high spatial/spectral res¬ 
olutions demonstrate that spiral galaxies often exhibit 
smooth ve locity gradients (Daigle et al. 2006 Dicaire 
et al.| 20 () 8 ), whereas mer gers tend to show m ore compli¬ 
cated kinem atic features (Mihos & Bothun||1998 Colina 
et al.||2005). Various visual and automatic analysis have 
also been developed to classify the me rger/disk nature of 
1 — 3 star-forming galaxies (e.g. Flores et al. 2006[ 


Shapiro et al.||2008[ [Hammer et al.||20b9[ [Bellocchi et al'l 

2012 Swinbank et al. 2012|. 


Large and systematic studies of the kinematic proper¬ 
ties of galaxies have been made possible with recent IFS 
surveys of nea rby and distant ga l axies (e.g. CALIFA sur¬ 
vey at z ^ 0: Husemann et a l. (2013); IMAGES survey 
at z ^ 0.6: [F l ores et al.|([^ 06|); MASSlV survey at z ^ 1: 
Epinat et ah (|200gt )1012|), SINS survey at z ~ 2: Forster 


Schreiber et al. 12006’3009)). The samples of galaxies 


with resolved kinematic information will be further ex¬ 
panded with the ongoing and future IFS sury e ys such as 


2014 


Wisnioski 
— 3 star- 


MaNGAFI a nd KMOS^^ (jWuyts et ah, 
et al.||20M|). Recent IFS surveys toward z 
tormmg galaxies (with typical star formation rates (SFR) 
of ^ 10-500 Mq yr“^) find that ~ 30 — 50% of these 
galaxies show evidence of interactions, yet a significant 
population of these star-forming galaxies (~ 40 — 70%) 


have kinematic pro 

parties consistent with a rotating disk 

(e.g. [Shapiro et al. 

L-iMi.s 

Law et al. 

2009 

Wright et al. 

2009[ Epinat et al. 

2012). 


A^ossible explanation of the large fraction of rotat- 


through rapid, smooth gas accretion (e.g. 

Keres et al. 

2005 Birnboim et al. 12007 

Dekel et al. 2009b). However, 


reformed in the merger remnant (Barnesj 

2002 

Sprinpel 

& Hernquist||2005l [Robertson & BullockI 

W 

Hopkins 

et al.||2009 Narayanan et al. |20(J9p. in tact, interacting 


tures along the merge r sequence (Mihos fc Bothun||1998 


Bellocchi et al. 2013), and it is unclear what the like¬ 
lihood is that mergers may be mistake nly classified as 
disk g alaxies or vice versa. For instance, [Bellocchi et al. 
( 2012 ) demonstrate that post-coalescence mergers can be 


miss-classified as disks whe n using the kinematic classi¬ 
fication criteria designed by Shapiro et al. (2008). 

In this paper, we further explore the robustness of var¬ 
ious kinematic classifications for interacting galaxies us¬ 
ing a sample of local (U)LIRGs spanning a wide range 
of merger stages and isolated disk systems. To make 
a realistic comparison with the observations of high—z 
galaxies, we artificially redshift our local sample as if 
they are observed at z = 1.5. We describe our sample 
and dataset in Section 2. In Section 3, we describe our 
analysis including the redshifitng procedures, the spec¬ 
tral line fitting, and the kinematic classification schemes. 
We present our results in Section 4 and discuss their im¬ 
plications in Section 5. We list our conclusions in Section 
6 . Throughout this paper, we adopt a ACDM cosmology 
with jJn = 70 k m s“^ Mpc“^, Dm = 0.3 and Da = 0.7 
( Hinshaw et al.|[2009 ). 


2. SAMPLE AND DATA 
2 . 1 . IFS Datacubes 

The IFS data of the 24 non-interacting and interacting 
systems used in this paper are obtained from the Wide 
Field Spectrosgraph (WiFeS), Integral Field Unit (IFU), 
and Great Obser vatory Allsky LIRGs Survey (G OALS) 


sample (WIGS: |Rich et al.j 2011 


GOALS survey (jArmus et al. 2009) comprises 203 ul- 


[2M4| . The 


traluminous and luminous infrared galaxies ((U)LIRGs; 
Air> lO^^L 0 )[jat z < 0.088, which is a complete sub- 
set of the IRAS Rev ised Bright Galaxy Sample (RBGS: 
Sanders et al.||2003|). The WIGS is a subset of GOALS 
galaxies that covers a wide range of merger stage s and 
luminositie s with a declination upper limit of -1-15° (Rich 


Lhe o bservations and data reductio n are described in 
detail in Rich et al.j (2010 [2011 2012), and here we pro¬ 
vide a brief summary. ' The IFS data were taken in 2009 
and 2010 using WiFeS ( ]Dopita et al.||2007 2010) on the 
ANU 2.3-meter telescope at the Siding Spring Observa- 
tory, and the typical seeing during the observations was 
^ 1.5". Our analysis is restricted to the red channel, 
with datacubes taken at a spectral resolution of R=7000 
and a wavelength coverage of ~ 5500 — 7000 A. Each 
WiFeS observation consists of 25 slitlets (l"x 38"). The 
WiFeS detectors have a pixel scale of 0.5" along the slit, 
and the final datacubes are binned by a factor of 2 so 
that the final pixel size is 1" x 1" (~ 400 pc at the me¬ 
dian z = 0.02). The typical total exposure time for each 
galaxy is 20—30 minutes. All of the final datacubes in 
this sample consist of at least two individual observations 
and in some cases a mosaic of multiple pointings to cover 
the extended structures. 

2 . 2 . Morphological Classifications 

The interaction stages of the 24 WIGS sources (non¬ 
interacting galaxies and interacting systems) were deter¬ 
mined by inspection of their optical morphology (Larson 
et al. 2014, in prep.). The visual classification scheme 
developed by Larson et al. (2014) is based on the mor¬ 
phological features d uring the progression of a n inter¬ 
action sequence (e.g. [Barnes & Hernquist||1992D: single 
disk galaxies (S), galaxy pairs without signs of interac- 
tions (Ml), galaxy pairs with signs of interactions (M2), 
merged galaxies with double nuclei (M3), merged galax¬ 
ies with single nucleus and tidal tails (M4), and merged 
galaxies with an offset single nucleus but without obvious 
tidal tails (M5). Go mparing to the visual classi fication 
scheme dev eloped by Veilleux et al. (2002) and Surace 
et al. (1998), Ml stage includes botn First ap proach and 
First contact s tages in Veilleux et al. (2002). The Pre¬ 
merger stage in Veilleux et ai. ] 2002 ) is divided into M2 or 
M3 depending on whether the disk-structure of galaxies 
can still be seen. The revised scheme allows us to probe 
the galaxy kinematics before and after galaxy disks being 
destroyed during the interaction. 

Eighteen of the WIGS sources have been observed with 
HST ACS in the E435kF band and FilAW band with 


^ https://www.sdss3.org/future/manga.php 


a resolution of 0".05 (Kim et al. 2013 Evans et al. in 
prep.), where the high resolution and sensitive HST im¬ 
ages can resolve the detailed structures and reveal the 

^ Ljr = Ls— iooOaht* ™ the object’s rest-frame 
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TABLE 1 

List of WIGS Sources 


IRAS Name 
(1) 

RA (J2000) 
(2) 

Dec (J2000) 
(3) 

Z 

(4) 

log(I/IR/L0) 

(5) 

Image/Filter 

(6) 

Int. Stage 

(7) 

Other Names 
(8) 

F01053-1746 

01:07:47.18 

-17:30:25.3 

0.020067 

11.71 

ACS F814W 

M3 

IC 1623, Arp 236 

F06076-2139 

06:09:45.81 

-21:40:23.7 

0.037446 

11.65 

ACS F814W 

M2 


08355-4944 

08:37:01.82 

-49:54:30.2 

0.025898 

11.62 

ACS F814W 

M3 


F10038-3338 

10:06:04.80 

-33:53:15.0 

0.034100 

11.78 

ACS F814W 

M3 


F10257-4339 

10:27:51.27 

-43:54:13.8 

0.009354 

11.64 

ACS F814W 

M4 

NGC 3256 

F12043-3140 

12:06:51.94 

-31:56:58.5 

0.023203 

11.36 

DSS R 

M2 


F12592-I-0436 

13:01:50.80 

-1-04:20:00.0 

0.037483 

11.68 

ACS F814W 

M4 


13120-5453 

13:15:06.35 

-55:09:22.7 

0.030761 

12.32 

ACS F814W 

M4 


F13373-I-0105 

13:39:55.00 

-1-00:50:07.0 

0.022500 

11.62 

ACS F814W 

M2 

NGC 5257/5258, Arp 240 

F15107-I-0724 

15:13:13.09 

-1-07:13:31.8 

0.012999 

11.35 

UH-2.2m I 

S 


F16164-0746 

16:19:11.79 

-07:54:02.8 

0.027152 

11.62 

ACS F814W 

M4 


F16399-0937 

16:42:40.21 

-09:43:14.4 

0.027012 

11.63 

ACS F814W 

M3 


F16443-2915 

16:47:31.06 

-29:21:21.6 

0.020881 

11.37 

IRAC 3.6 fim 

S 


F17138-1017 

17:16:35.79 

-10:20:39.4 

0.017335 

11.49 

ACS F814W 

S 


F17207-0014 

17:23:21.95 

-00:17:00.9 

0.042810 

12.46 

ACS F814W 

M4 


F17222-5953 

17:26:43.34 

-59:55:55.3 

0.020781 

11.41 

DSS R 

S 


F18093-5744 

18:13:39.63 

-57:43:31.3 

0.017345 

11.62 

ACS F814W 

M2 

IC 4686/4687/4689 

F18293-3413 

18:32:41.13 

-34:11:27.5 

0.018176 

11.88 

ACS F814W 

M2 


F18341-5732 

18:38:25.70 

-57:29:25.6 

0.015611 

11.35 

DSS R 

S 

IC 4734 

F19115-2124 

19:14:30.90 

-21:19:07.0 

0.048727 

11.93 

ACS F814W 

M2 


F20551-4250 

20:58:26.79 

-42:39:00.3 

0.042996 

12.06 

ACS F814W 

M4 


F21330-3846 

21:36:10.83 

-38:32:37.9 

0.019060 

11.14 

DSS R 

M2 


F22467-4906 

22:49:39.87 

-48:50:58.1 

0.043033 

11.84 

ACS F814W 

M4 


F23128-5919 

23:15:46.78 

-59:03:15.6 

0.044601 

12.06 

ACS F814W 

M3 



NOTE: Columns (1) — (5) list the IRAS name, RA, Dec, 2 , and for the WIGS sources, which are compiled from|Sanders et al.|||2003|l 
and NASA/IPAC Extragalactic Database (NED). Column (6) indicates the images and filters that are used for visual morphological 
classification, column (7) lists the morphological classes from Larson et al. (2014, in prep.), and column (8) lists other common names of 
WIGS sources. 


faint interacting features. The F814tT band images are 
used for morphological classification due to their less se¬ 
vere dust obscuration than the F4351T band images. In 
fact, dust obscuration may remain an issue at FilAW 
band, which can bias our classihcation at M3-M4 stages. 
Nine out of 12 sources that are classified as M3 or M4 
in our sa mple have NICMOS F160W band images avail¬ 
able from Haan et al. (2011). The classification results 
of these 9 sour ces do not clmnge when f 1601T ban d im¬ 
ages are used (|Haan et al.||2011| |Petty et al.||2014|). For 
the other 6 sources that do not have JUs T linages avail¬ 
able, we use R band images from Digitized Sky Survey 
(DSS), I ban d images from the UH-2.2 meter telescope 
(Ishida 2004), or post calibration Spitzer IKAC 3.6pm 
images from NASA/IPAC Infrared Science Archive to 
determine their morphological classes. The plate scales 
of DSS, UH 2.2-meter, IRAC images used in this paper 
are 1", 0".22 (seeing-0".9), and 0".09 (FWHM~1".66), 
respectively. Based on these images, we have 6, 6, 5, 7 
sources classihed as S, M2, M3, M4, respectively. Ta¬ 
ble 1^ lists basic information of galaxies, optical images 
usea for morphological classifications, and the interac¬ 
tion stages. We note that three WIGS sources include 
multiple components that are individually resolved in the 
IFS observations, and we carry out their kinematic clas¬ 
sifications separately. The kinematic classifications of a 
total of 28 galaxies are listed in Table and their optical 
images are shown in Figure 

3. ANALYSIS 

3.1. Redshifted IFS Datacubes 


To test if the merger/disk classification based on kine¬ 
matic properties can robustly identify high—z interacting 
galaxies, we artificially redshift the IFS datacubes of lo- 


surveys cover a range 

of redshifts from z ^ 1 — 3 (Forster! 

Schreiber et al. 120091 

Law et al.|2009|| Wright et al.|2UUyt 

Wisnioski et al.||2014 

L h'or example, the SINk'UNI sam- 


1.2 (Epinat et al. 


2012). Recent KMOS^^ 


2009 

survey focuses on star-iorm ing galaxies at z ~ 1 and 
z ^ 2 ( Wisnioski et al.|2014 ), and the AMAZE and LSD 


3 ( IGnerucci et al.pOll ). 
edsnitt our local galaxies 


proj ects target galaxies at z 

In this paper, we artificially redshift our local galaxies 
to z = 1.5, as a case study. The redshift is chosen to 
match a lensed spiral galaxy at z = 1.5 (Spll49), whose 
fully reduced Keck/OSIRIS observations are available to 
simulate the sky- line dominated n oise of near infrared 
IFS observations (Yuan et al.||2011 ). The corresponding 
physical sizes of a given angular scale only differ by ~ 1% 
between z = 1.5 and z = 2. Thus our choice of z = 1.5 
should be representative of typical z ^ 1 — 2 observations 
in terms of the spatial resolutions can be achieved. De¬ 
tail investigations of the redshift dependence regarding 
to the sky line distributions and instrument conhgura- 
tions are beyond the scope of this paper. We assume the 
observations are carried out under seeing-limited condi¬ 
tion (0".5) with pixel size of 0".l x 0".l 900 pc at 

z = 1.5). 

To redshift the local galaxies, we first define an el¬ 
liptical mask for each galaxy. This is to prevent con¬ 
tamination from foreground stars and other artifacts 




































4 


when binning the pixels and convolving with larger 
beams. We then follow the r edshifting procedu re of 
IRASF17222—5953 described in |Yuan et al.| (|2013[ Sec¬ 
tion 4.2). We derive the surface brightness distribution 
according to the redshift and plate scales of the input 
(WiFeS) and output (OSIRIS) datacubes. For instance, 
the WiFeS plate scale (l"xl") corresponds to ~400 pc 
X 400 pc at 0 = 0.02, and we then bin the pixels by a 
factor of two to match the corresponding physical sizes 
of the 0".l X 0".l plate scale at z = 1.5 (~900 pc x 900 
pc). 

Then we use the source-free region of the Spll49 data 
to create noise datacubes with the dimensions of the 
spatially binned and dimmed WiFeS data. The noise 
datacube is interpolated to match the spectral resolu¬ 
tion of WiFeS data. We combine the spatially-binned 
WiFeS data with the noise datacube, and then convolve 
the resulting datacube with a two-dimensional Gaussian 
(FWHM=0".5), a resolution comparable to the typical 
seeing-limited observations of z ~ 1 — 3 star-forming 
galaxies. Since the emission lines of these degraded data 
become too faint to be detectable, we apply an artificial 
brightening factor to ensure that the peak of Ha emission 
has a signal-to-noise (S/N) ratio of 50. 

3.2. Spectral Line Fitting 

To extract the kinematic information from the artifi¬ 
cially redshifted datacube, we fit Gaussian profiles to Ha 
{\rest-frame = 6562.8 A and two adjacent [N ii] lines 
{Xrest-frame = 6548,6583 A simultaneously. We fix the 
flux ratio betw een [Nii]A6583 and [Nii]A6548 at 3 (|Os- 
terbrock||1989|) , and we assume the linewidths of the Ha 
and [IN ii] lines are the same. An instrumental profile of 
0.94 A is then subtracted from the measured linewidth 
and the final linewidth = “ ^Lstmmentf'^ ■ 

In the end, our fitting profile contains five free param¬ 
eters: systematic velocity, spectral linewidth, Ha flux, 
[N ii]A6548 flux, and the continuum level (A6400 — 6750 

A). 

We use t he IDL routine MP FIT to perform the spectral 
line fitting ( Markwardt|2009 ), which uses the Levenberg- 
Marquardt afgorithm to solve non-linear least-squares 
problems. We repeat the spectral line htting at each pixel 
across the entire map, and each line htting is weighted 
based on the inverse squared noise spectrum. The un¬ 
certainties of each free parameter are the 1 a errors com¬ 
puted from the covariance matrix in the MPFIT routine. 
We use the Ha line as a tracer of a galaxy’s kinemat¬ 
ics. A minimum of S/N=5 is required for the Ha line t o 
be considered a robust detection (e.g. Yuan et n!||2013 ). 
No signihcant differences are seen in the bins when us - 
ing adaptive spatial binning ( Cappellari fc Copin||2003 ). 
Based on the best ht systematic velocity and line width 
at each pixel, we then construct the velocity distribu¬ 
tion map and the velocity dispersion map. An example 
galaxy (IRASF01053—1746) is shown in Figure]^ and 
we present the entire sample in the Appendix (Figure]^. 

3.3. Kinematic Classifications 

Several kinematic classihcation schemes have been de¬ 
veloped to classify merger/no n-merger of local and z 
1—3 star-forming galaxies (e.g. [Flores et al.|2006[ [Shapiro 


20121 [Bellocchi et al.|[20T2| |Epinat et al.|[20T2{ |Belloc-| 
chi et al. 2013p. Here we use the classihcation schemes 


based on ki nematic asymmetries develo ped by Shu|)jjo 


et al. 


(|2008[ hereafter S08) and revised by Bellocchi et al. 


)2U12[ hereafterB12), as well as a modihed version of the 


visual classihcation scheme developed by |Flores et al.| 

(l2006|. 


3.3.1. Kinematic Asymmetries in\Shapiro et al.\ l[200S\) 


The kinemetry analysis (Krajnovic et al. 2006) is de¬ 
signed for modeling the higher-order moments of the ve¬ 
locity and velocity dispersion distributions of galaxies. 
The line-of-sight velocity map or velocity dispersion map 
K{a,ijj) can be divided into several elliptical rings (with 
semi-major axis a) as velocity or velocity dispersion pro- 
hles. These prohles can then be described as an expan¬ 
sion of N + 1 harmonic terms: 


N 


K{a,il)) = Ao(a)-|- > A„(a) sinn'!/'-|-H„(a) cosnr/>, (1) 


where '0 is the azimuthal angle. For an ideal disk, its ve¬ 
locity held and velocity dispersion held are dominated by 
Bi and Aq, and thus higher-order deviations may be a re¬ 
sult of galaxy mergers. S08 quantify the level of deviation 
from an ideal disk by dehning asymmetric measures of 
velocity and velocity dispersion helds [vasym and <Jasym) 
using the higher order kinematic coefficients (Equation 4 
and 5 in SOS). Based on the templates from local galaxies 
and simulated galaxies, SOS dehne a criterion to identify 
mergers via Ka — _l .^2 


^asym (^asym ^asym) ^ ^ 0.5. 


We measure Vasym and aasym of our redshifted WIGS 


et al.|200^|Gonqarves et al.|2l)I0[|Alaghband-Zadeh et al. 


galaxies us ing the IDL routine Kinemetry (Krajnovic 
et al. 2006). Before performing the kinemetry analysis, 
we determine the position of the galaxy center based on 
the distribution of the continuum emission. The contin¬ 
uum map is hrst smoothed using a Gaussian with FWHM 
of 3 pixels (corresponds to 0".3 in the redshifted dat¬ 
acubes), and then the galaxy center is determined as the 
centroid of the smoothed map (typical difference of us¬ 
ing different smoothing factors of 1-5 pixels is less than 
1 pixel). We note that the smoothed continuum maps 
are only used for determining galaxy center. We avoid 
using the Ha emission since it traces the recent star form¬ 
ing regions and may not necessarily represent the overall 
stellar distribution. As discussed in S08 and B12, the 
derivation of kinematic asymmetries is sensitive to the 
position of the galaxy center determination. We discuss 
how the accuracy of the galaxy center impacts our results 
and conclusions in Section 4.3. 

After fixing the center, kinemetry then finds the best- 
fit ellipse with position angle (PA) and the flattening fac¬ 
tor (Q = 1 —e) at each radius step until more than 25% of 
the data points along an ellipse are not present. S08 de¬ 
rive a global PA and Q before performing the kinemetry 
analysis, but here we treat them as free parameters since 
the meaning of a global PA and Q for mergers is unclear. 
B12 also show that the choice of either free or fixed PA 
or Q does not have a strong impact on their results. We 
derive Kasym of six single disk galaxies by attempting 
to find a global PA and Q prior to the kinemetry anal¬ 
ysis, none but one galaxy (F16443—2915 S) alters the 
classification result. 
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Fig. 1.— The left two panels shows the ACS/F’814VF-band image of IRASF01053—1746 with two different stretching schemes (inverse 
hyperbolic sine and histogram equalization) and its morphological classification (M3). The right two panels show the velocity map and 
velocity dispersion map derived from the redshifted IFS datacubes. The velocity and velocity dispersion maps have the same field-of-view 
in physical scale as the optical images shown in the left panel, and their angular sizes are scaled to z = 1.5. The color bars are in the unit 
of km s“^. We present the maps for the entire sample in the Appendix (Figure]^. 


TABLE 2 

Results of Kinematic Classifications 


IRAS Name 

Int. Stage 

Kinemetry in SOS 
I^asym (Classification) 

Revised Kinemetry in B12 
Kasym. B12 (Classification) 

Visual Kinematic 
Classification 

E01053-1746 

M3 

0.16±0.04 (Disk) 

0.20±0.05 (Merger) 

RD 

E06076-2139 

M2 

1.82±0.52 (Merger) 

1.38±0.43 (Merger) 

CK 

08355-4944 

M3 

4.06ibl.80 (Merger) 

4.39±1.81 Aerger) 

CK 

F10038-3338 

M3 

5.18±2.39 (Merger) 

7.23±2.92 (Merger) 

CK 

F10257-4339 

M4 

0.05±0.01 (Disk) 

0.04±0.03 (Disk) 

RD 

F12043-3140 

M2 

0.27±0.07 (Disk) 

0.34ib0.07 (Merger) 

PD 

E12592-I-0436 

M4 

0.15±0.01 (Disk) 

0.12±0.01 (Disk) 

RD 

13120-5453 

M4 

0.39±0.08 (Disk) 

0.33±0.07 (Merger) 

RD 

E13373-I-0105 W 

M2 

0.09±0.01 (Disk) 

0.07±0.01 (Disk) 

RD 

E13373-I-0105 E 

M2 

0.14±0.02 (Disk) 

0.12±0.01 (Disk) 

RD 

F15107-I-0724 

S 

0.50±0.06 (Disk) 

0.52±0.05 (Merger) 

PD 

F16164-0746 

M4 

2.31±0.77 (Merger) 

2.38±0.77 (Merger) 

CK 

F16399-0937 

M3 

0.21±0.06 (Disk) 

0.20ib0.06 (Merger) 

PD 

F16443-2915 S 

S 

0.34±0.06 (Disk) 

0.41di0.07 (Merger) 

PD 

E16443-2915 N 

S 

O.lOibO.Ol (Disk) 

0.07±0.01 (Disk) 

RD 

E17138-1017 

s 

0.29±0.04 (Disk) 

0.27±0.03 (Merger) 

PD 

E17207-0014 

M4 

0.19±0.08 (Disk) 

0.26±0.10 (Merger) 

PD 

E17222-5953 

S 

0.11±0.02 (Disk) 

0.13±0.02 (Disk) 

RD 

F18093-5744 N 

M2 

0.11±0.02 (Disk) 

O.llibO.Ol (Disk) 

RD 

F18093-5744 S 

M2 

0.18±0.03 (Disk) 

0.21ib0.04 (Merger) 

RD 

F18093-5744 C 

M2 

0.76±0.10 (Merger) 

0.55±0.12 (Merger) 

CK 

F18293-3413 

M2 

0.24±0.06 (Disk) 

0.26±0.07 (Merger) 

PD 

E18341-5732 

S 

O.lOibO.Ol (Disk) 

0.09±0.01 (Disk) 

RD 

E19115-2124 

M2 

0.57±0.24 (Merger) 

0.93± 0.32 (Merger) 

CK 

E20551-4250 

M4 

0.31±0.03 (Disk) 

0.25±0.02 (Merger) 

PD 

F21330-3846 

M2 

0.41±0.24 (Disk) 

0.60ib0.31 (Merger) 

PD 

F22467-4906 

M4 

0.82±0.36 (Merger) 

0.62±0.27 (Merger) 

CK 

F23128-5919 

M3 

0.43±0.06 (Disk) 

0.48±0.07 (Merger) 

CK 


NOTE: Kinematic classifications of S, M2. M3, and M4 galax ies based on the classification schemes described in Section 3.3.1- 3. The 
merger/disk classification criteria in the|Shapiro et al.|(|2008j scheme follows Kasym > 0.5, and the classification criteria in the|Bellocchi | 
|et aL| |2012| follows Kasym, B\2 > 0.146. 


3.3.2. Revised Kinematic Asymmetries in\Bellocchi et al.\ 


Although the local merger templates used in SOS 
include 3 interacting pairs, 2 double nuclei mergers, 
and 3 single nucleus mergers, the sample size remains 
small and thus they may not trace the full scheme 
of the interaction sequence. B12 find that their sam¬ 
ple of both non-interacting disks (IRASF10567—4310 
and IRASF11255—4120) and post-coalescence mergers 


(IRASF04315-0840 and IRASF21453-3511) are classi¬ 
fied as disks using the kinematic classification scheme 
developed by SOS. Since a post-coalescence merger may 
be dominated by rotation in the inner regions but re¬ 


tain distorted kine matics at the outskirts (e.g. Kron- 


berger et al. 2007), B12 modify the asymmetric mea- 
sures defined in SOS by weighting the asymmetries ac¬ 
cording to the circumferences of each radius step (Equa¬ 
tion (6) and (7) in B12; hereafter indicated as Vasym,Bi 2 
and <Jasym,Bi 2 )- Based on the revised measurements of 
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kinematic asymmetries of two non-interacting disks and 
two post-coalescence mergers, B12 derive a new crite¬ 
rion to classify mergers via Kasym.,B 12 = iylsym,Bi 2 + 

<^Lym,B12y^‘^ > 0 - 146 . 


3.3.3. Visual Kinematic Classification 

In addition to the two automatic classification schemes 
defined by SOS and B12, we also visually classify the dy¬ 
namical status of galaxies based on the distribution of 
their velocity and velocity dispersion fields. We adopt 
the 3-class k i nemat ic classi fication scheme develop ed by 
Flores et al. (2006 also see Bellocchi et a l. (20 13t) with 
some modific ations. InjFlores et ai. 1 2006) and" Bellocchi 
et al. (2013), comparison between kinematic maps and 
optical morphology (alignments between optical major 
axis and velocity gradient) is required to determine the 
kinematic classes of galaxies. To make a direct compar¬ 
ison with the kinemetry-based classifications, which do 
not require the alignment with optical m ajor axis, we 
modify the visual classification s chemes in Flores et al.| 
(2006) and Bellocchi et al. (|2013 ) by removing the crite¬ 
ria that use optical morphology. 

Below we list the modified classification criteria of 
three kinematic classes determined based on a qualita¬ 
tive assessment of ralaxies’ velocity and velocity disper¬ 
sion maps (Figure!^: 


1. Rotating disk (RD): the velocity map shows a clear 
velocity gradient, and the velocity dispersion map 
should show a clear peak near the galaxy center 
where the gr adient of the rotation curve is the 
steepest (e.g. van Zee fc Bryant]|1999|. 


2. Perturbed disk (PD): the velocity map shows a 
clear velocity gradient as RD whereas the velocity 
dispersion map shows a peak offset from the center 
or does not show a clear peak. 


3. Complex kinematics (CK): both velocity and ve¬ 
locity dispersion maps do not follow the expected 
pattern of normal rotating disks. 


4. RESULTS 

4.1. Kinemetry-hased classifications along the 
interaction sequence 

All the kinematic classifications of S, M2, M3, and 
M4 galaxies based on the classification schemes in Sec¬ 
tion 3.3 are listed in Table Figure shows the dis¬ 
tribution of kinematic asymmetries denned by SOS and 
B12 for our WIGS galaxies, which are artificially red- 
shifted to z = 1.5 assuming seeing-limited conditions 
(FWHM=0".5). The data points are color-coded accord¬ 
ing to their interaction stages and the error bars rep¬ 
resent the standard deviation of the mean Vasym and 
c^asym across all radius steps. When the SOS scheme 
is used, all 6 single disk galaxies are classified as disks 
based on the kinematic asymmetries, yet only S out of 
22 interacting systems are classified as mergers. As an 
example, IFASF12592-I-0436 demonstrates the case that 
the morphological and kinematic classifications are in¬ 
consistent. Despite the obvious tidal tails in the F814IF- 
band image, the velocity and velocity dispersion fields of 


IFASF12592-I-0436 do not show complicated structures 
as one might expect for a merger. 

The consistency between merger/disk classifications 
based on morphological features and kinematic asymme¬ 
tries shows a strong trend with galaxy interaction stage. 
Figure shows the fraction of galaxies that are classified 
as disks and mergers based on kinematic asymmetries 
as a function of morphological classifications. When the 
classification criteria of SOS is used, 100% of single disk 
galaxies in our sample are classified as disks, whereas 
only 30-40% of the interacting systems are classified as 
mergers. Since the interacting disks (M2) still show clear 
disk structures in spite of their signs of interactions, these 
galaxies may retain disk-like kinematics. However, our 
sample of later-stage mergers (mergers with double nuclei 
and single nucleus; M3 and M4) display highly disturbed 
morphological features, yet only ~ 30 — 40% are classified 
as mergers. 

On the other hand, when the classification scheme 
based on the revised kinematic asymmetries in B12 is 
used, all of our M3 galaxies are classified as mergers, 
whereas the success rate for single disk galaxies to be 
classified as disks decreases to only ~ 50%. The dif¬ 
ferences in the definition of kinematic asymmetries (by 
simply averaging across all radius steps or weighted ac¬ 
cording to the circumferences) may contribute to the dif¬ 
ferences of merger/disk fractions measured using the SOS 
or B12 scheme. However, the median difference of the 
Kasym and Kasym,B 12 values is only ^ 16% for our galax¬ 
ies (Table [^, whereas the classification criteria used in 
B12 (ATasym,512=0.146) is more than three times lower 
compared to SOS {Kasym=0-5). Thus the significantly 
lower classification threshold used in B12 is the main 
reason of the drastically different disk and merger frac¬ 
tions derived using the SOS and B12 schemes. We note 
that B12 shows that the classification criterion between 
mergers and disks decreases from 0.146 to 0.13 when ar¬ 
tificially redshifting their datacubes to z = 3, and thus 
a lower threshold than 0.146 should be applied to our 
redshifted datacubes. However, none of our classifica¬ 
tion results changes using a criterion of either 0.146 or 
0.13. Figure demonstrates that there is a significant 
overlap between interacting systems and single galaxies 
in the distribution of kinematic asymmetries using either 
SOS or B12 scheme, which implies that a single Kasym or 
Kasym,B12 cut is Unable to differentiate interacting sys¬ 
tems from single galaxies. A higher threshold can lead to 
an underestimation of merger fractions whereas a lower 
threshold can overestimate merger fractions. 

The results on the reliability of merger/disk classifica¬ 
tions based on kinematic properties rely on robust mor¬ 
phological classifications. Six galaxies (4 S galaxies, 2 M2 
galaxies) included for kinematic classifications have only 
low resolution DSS or IRAC images available (with the 
resolution >1"), and their morphological classifications 
are less reliable compared to the other 22 galaxies. If we 
exclude the six galaxies with low resolution images, only 
2 S galaxies remain in our sample. Both of these two S 
galaxies are classified as disks when the SOS scheme is 
used, and are classified mergers when the B12 scheme 
is sued. For M2 galaxies, the merger fraction increases 
from 30% to 37.5% when the SOS scheme is used, and 
decreases from 70% to 62.5% when the B12 scheme is 
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®asym,B12 


Fig. 2.— The distribution of kinematic asymmetries for our redshifted WIGS galaxies (with a spatial resolution of 0”.5). The left panel 
shows the results based on Vasym and cFasym defined by SOS. Data points are overlaid on the merger (red) and disk (blue) templates 
distribution in SOS, and the black dotted line refers to the merger/disk classification criteria of Kasym = (i^asym + ^asym)^^^ — 0-5. 
The right panel shows the results based on the revised asymmetries defined by B12, and the black dotted line refers to the merger/disk 
classification criteria of Kasym,B 12 = ('’^asym B 12 ^asym 312 )^^^ “ 0.146. In both panels, data points are color-coded according to their 
morphological classifications, in which the classification scheme is defined based on the progression of an interaction sequence: gray (S: 
single disk galaxies), cyan (M2: interacting disks), orange (M3: mergers with double nuclei and tidal tails), and yellow (M4: mergers with 
single nucleus and tidal tails). The error bars represent the standard deviation of the mean Vasym and Oasym across all radius steps. 
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Fig. 3.— The fraction of galaxies to be classified as disks/mergers 
based on the measurement of kinematic asymmetries as a function 
of interaction stages (S: single disk galaxies, M2: interacting disks, 
M3: mergers with double nuclei, and M4: mergers with single nu¬ 
cleus) . The black dots are the results based on the criteria defined 
by S08 and the red dots are based on the criteria defined by B12. 
Error bars are derived based on the frequency (out of four addi¬ 
tional testing cases) of altering the merger/disk classification when 
shifting the center of galaxies horizontally and vertically by 1 pixel. 


used. 


4 . 2 . Comparison between visual kinematics and 
kinemetry-based classifications 

Figure 1^ shows the classification results based on the 
visual kinematic classification described in Section 3.3.3. 
The distribution of rotating disks (RD), perturbed disks 
(PD), and complex kinematics (CK) along the inter¬ 
action stages is consistent with the trend based on 
kinemetry-based classifications. All of our single disk 
galaxies are classified as either RD or PD. As the in¬ 
teraction stage progresses, higher fraction of the galax¬ 


ies begin to show complicated kinematic features. For 
instance, 60% of the M3 galaxies are classified as CK. 
However, a significant population of M2, M3, and M4 
galaxies retain an ordered velocity gradient in their ve¬ 
locity maps and thus are classified as either RD or PD. 
We compare our visual k inematic classification results 
with Bellocchi et al. (2013) based on the overlapped sam¬ 
ple of 11 galaxies. As noted in Section 3.3.3, the vi sual 
classification scheme used in Bellocchi et al. (2013) re¬ 
quires the comparison with optical major axis. In spite 
of this difference and that our redshifted datacubes have 
a factor of 2- 6 worse spatial resoluti on in physical scales 
compared to Bellocchi et al. (2013), seven out of these 
11 galaxies have consistent visual kinematic classifica¬ 
tion results, and the discrepancies in the other 4 galaxies 
are mainly due to the degree or complexity (PD or CK) 
or the requirement of alignment with optical maior axis 
(RD or PD). 

Table [^compares the results between visual kinematics 
and kinemetry-based classifications. When SOS criteria 
is used, all galaxies identified as RD are classified as disks 
and all but one galaxies identified as CK are classified as 
mergers. When B12 criteria is used, all but three identi¬ 
fied as RD are classified as disks and all galaxies identi¬ 
fied as CK are classified as mergers. The major difference 
is those galaxies that are identified as PD, where all of 
them are classified as disks in SOS criteria but classified 
as mergers in B12 criteria. The comparison results show 
that different kinematic classification schemes can lead to 
large dispersions in the population of identified mergers. 


4 . 3 . Uncertainties in the kinemetry-based classifications 

As mentioned in Section 3.3.1, the derivation of Vasym 
and CFasym IS Sensitive to the position of the galaxy center. 
The uncertainties of galaxy center can be a significant 
factor when different determination methods are used. 
For instance, whether we use the position of the peak 
































5. DISCUSSION 



Interaction Stage 

Fig. 4.— The fraction of S, M2, M3, and M4 galaxies to be clas¬ 
sified as rotating disks (RD), perturbed disks (PD), and complex 
kinematics (CK). 


TABLE 3 

Comparison between visual kinematics and kinemetry-based 
CLASSIFICATIONS 


Visual Classes 

RD 

PD 

CK 

SOS Disks 

11 

9 

1 

S08 Mergers 

0 

0 

7 

B12 Disks 

8 

0 

0 

B12 Mergers 

3 

9 

8 


emission or derive the centroid with different smoothing 
factors (see the description in Section 3.3.1). We derive 
the difference between the position of the center deter¬ 
mined in Section 3.3.1 and the cases that are determined 
based on the peak of continuum maps or another four 
different smoothing factors (1,2,4,5 pixels), and we find 
that the average difference is 0.45 pixel (with the stan¬ 
dard deviation of 0.44 pixel). In fact, one of our WIGS 
sample overlaps with the merger templates used in SOS 
(IRASF17207—0014), yet our derived Kasym is ^ 6 times 
smaller compared to the median Kasym of the Monte 
Carlo realizations in SOS. Although the differences be¬ 
tween the datacubes used in SOS and in this paper may 
explain the difference in Kasym, the most likely expla¬ 
nation is the difference in galaxy center (the center de¬ 
termined by SOS is located close to the edge of the IFU 
datacube; Figure 9 in SOS). A similar value of Kasym 
as SOS is derived when shifting the galaxy center by 1-2 
pixels. 

Since the pixels are discretized (each pixel corresponds 
to 0".l and ^0.9 kpc in the redshifted datacubes), we 
estimate the significance of galaxy center determination 
to our conclusions by calculating the kinematic asymme¬ 
tries in four additional cases for all galaxies assuming the 
galaxy center is shifted by one pixel horizontally or verti¬ 
cally with respect to the center determined based on the 
continuum map. We determine the uncertainties in frac¬ 
tion of galaxy classification results (Figure based on 
the likelihood that the merger/disk classification alters 
when shifting the galaxy center. The resulting uncer¬ 
tainties are typically 10 — 20%, and thus our conclusions 
remain after considering the uncertainty in the determi¬ 
nation of galaxy center. 


By carrying out kinematic classifications of a set of 24 
interacting systems and non-interacting galaxies, we find 
that mergers with two separate nuclei (M3 galaxies) are 
the most easily distinguished from single disk galaxies, 
and a higher fraction of early and late stage mergers (M2 
and M4) have kinematic properties that can be confused 
with rotating disks. 

The high fraction of interacting galaxies showing kine¬ 
matics consistent with disks is in disagreement with the 
conclusions based on local merger templates in SOS, and 
this differences cannot be explained by the uncertainties 
in the galaxy center determination alone. However, one 
local merger template (IRAS12112-I-0305, out of eight 
galaxies) in SOS that shows disk-like kinematics has been 
excluded for deriving the merger/disk classification cri¬ 
teria. If IRAS12112-I-0305 was included by SOS in de¬ 
termining merger/disk Kasym thresholds, the differences 
of the derived disk fractions between this work and SOS 
may be smaller. Furthermore, both of the WIGS and 
SOS samples are of modest sizes, and these samples may 
cover a different range of interaction sequence. For in¬ 
stance, we also include interacting pairs that are individ¬ 
ually resolved in the IFS observations, whereas the three 
interacting pairs used SOS are all unresolved. 

As discussed in Section 4.1, our M2 galaxies show clear 
disk structure in the optical images, and thus we ex¬ 
pect them to retain disk-like kinematics. Five out of 
six M2 systems that are individually resolved in the 
IFS observations F13373-b0105 W/E, FIS093-5744 N/S, 
FIS293—3413, exception: F18093—5744 G) show sig¬ 
nificantly lower kinematic asymmetries calculated us¬ 
ing either SOS or B12 schemes compared to the other 
four M2 systems that are unresolved in the IFS ob¬ 
servations (F06076-2139, F12043-3140, F19115-2124, 
F21330—3S46). This difference suggests that individu¬ 
ally resolved galaxies in the interacting pairs tend to 
show disk-like kinematics, whereas if the two galaxies 
are unresolved, then they tend to show more complicated 
kinematics. For example, the unresolved dynamical in¬ 
formation from two galaxies in F19115—2124 leads to its 
highly disturbed and complicated velocity field. 

One possible origin of the disk-like kinematics in the 


merger remnants (e.g. 

Barnes 

2002 

Springel & Hernquist 

2005 

Robertson & E 

lullock 

200S 

. Some of the local 

tujTi 

tiGs may have comparable gas tractions to those of 


that are necessary to reform the gas disk in galaxy merg - 
ers (e.g. |Downes & Solomon||199S| |Tacconi et al.||1999|). 
The emergence of these gas disks and their subsequent 
formation of stellar disks are also o bserved in luminous 
mergers and merger remn ants (e.g. Rothberg fc Fischer! 


2010[ Medling et al.|[20l4 |. Recently, Ueda et al. 1 20141 ) 
study the ^^GO kinematic pr operties of optically-s e lected 
local merger remnants from Rothberg & Joseph (20041 
(morphological properties consistent with our Md and 
M5 stages), and they find that S0% (24/30) of their sam¬ 
ple show kinematic signatures of rotating molecular gas 
disks. Such high disk fraction is consistent with the disk 
fraction in our M4 galaxies derived from the SOS classi¬ 
fication criteria and the visual classification schemes. 

Nonetheless, part of our M3 and M4 galaxies that show 
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disk-like kinematics may simply reflect the hete rogeneous 


dynamical phases during the interactions (e.g. Mihos & 

Bothun 1998 

Golina et al. |2005 iBellocchi et al.| 21)13). 

h'or instance. 

Bellocchi et al.|(|2L)13|) show that only about 


M3 and M4 in the definition of this paper) display com¬ 
plicated kinematics and the other half are classified as 
perturbed disks based on the visual classification scheme. 
In this case, higher spatial resolution observations have 
the potential to resolve detailed kinematic information 
and trace the disturbance due to galaxy mergers. Using 
the kinemetry-based classific ation scheme developed by 
SOS, iGongalves et al.| (|2010|) find that the merger frac¬ 
tion of a set of Lymaii Break Analogs (LBAs) at z ~ 0.2 
decreases from ^ 70% to ^ 38% after redshifting their 
sample to z = 2.2, which is partly due to the ~ 10 times 
worse resolution in the redshifted datasets. IFS obser¬ 
vations with high spatial resolution of local interacting 
systems will provide further insight on how spatial res¬ 
olution have impact on the kinemetry-based analysis at 
different merger stages. 

The significant overlap between single galaxies and in¬ 
teracting systems in the distribution of kinematic asym¬ 
metries (defined by simply averaging across all radius 
steps or weighted according to the circumferences) im¬ 
plies that a simple cut in Kasym or Kasym, bi 2 is unable 
to differentiate mergers from isolated disks. When the 
classification scheme in SOS is used, the success rate to 
classify S galaxies as disks is 100%, yet the success rate 
to classify M3 galaxies as mergers is only 40%. The B12 
scheme results in 100% success rate in classifying M3 
galaxies as mergers, whereas the success rate to classify 
S galaxies as disks decreases to only 50%. Meanwhile, 
we caution directly applying these merger recovery rate 
corrections to other work suc h as the samples in SOS or 


Alaghband-Zadeh et al. (2012) as this factor may vary de- 
pending on the physical properties of galaxy samples (e.g. 
merger stages, luminosity, mass, gas content). Other in¬ 
dicators such as multi-waveleng th morphological proper¬ 


ties (e.g. Hammer et al. 
constrain the merger/dis. 
forming galaxies. 


2009) are necessary to further 
c classification of z ^ 1 — 3 star- 


6 . CONCLUSIONS 

We have carried out merger/disk classifications for a 
set of 24 local (U)LIRGs based on the measurement of 
kinematic asymmetries and visual inspection. Our sam¬ 
ple spans a wide range of interaction sequence from sin¬ 
gle disk galaxies, interacting disks, to fully merged rem¬ 
nants. To make a realistic comparison with the observa¬ 
tions carried out for z ~ 1 — 3 star-forming galaxies, we 
degrade the spatial resolution and sensitivity of the IFS 
datacubes as if the galaxies are observed at z = 1.5 with 
a 0".5 seeing. The classifications of 28 local interact¬ 
ing systems and single disk galaxies lead to the following 


conclusions: 

1. Based on the merger/disk fractions derived from 
kinematic classifications, we find that mergers 
clearly displaying two nuclei and tidal tails (M3 
galaxies) show kinematic features most readily 
able to distinguish them from single disk galaxies. 
Nonetheless, some of our interacting disk pairs and 
merger remnants could be mistaken as disks when 
using kinemetry-based or visual classifications. 

2. We find that a high fraction of late-stage merg¬ 
ers (M4 galaxies) show signatures of rotating disks. 
This result may be explained as the disks reformed 
in the galaxy mergers or simply reflects the hetero¬ 
geneous dynamical phases during the interactions. 

3. Interacting disk pairs (M2 galaxies) that are re¬ 
solved in the IFS observations have, as expected, 
smaller values of both Kasym and Kasym, bi 2 com¬ 
pared to those pairs that are unresolved in the IFS 
observations. Therefore, during the early interac¬ 
tion stages, kinematic classifications are only sensi¬ 
tive to the systems with small projected distances. 

4. Our results show that kinematic properties alone 
are insufficient to constrain the merger fraction of 
z ^ 1 — 3 star-forming galaxies. Other merger indi¬ 
cators such as galaxy morphology traced by stars 
or molecular gas are required to further constrain 
the merger fraction at high—z. 
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APPENDIX 

A. OPTICAL IMAGES, Ha VELOCITY AND VELOCITY DISPERSION MAPS OF REDSHIFTED WIGS GALAXIES 

In this section, we present the optical images of the WIGS galaxies, velocity maps and velocity dispersion maps 
derived from the redshifted IFS datacubes. We include all 28 galaxies that are used for kinematic classifications. 
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High Contrast Image Low Contrast Image 

IRASF01053-1746 (M3) 


Velocity Map 


Velocity Dispersion Map 


30 “ 0 “ 




IRASF10257-4339 (M4) 

-100 0 100 30 60 90 120 



0 “ 0“.5 0 “ 0.5 


Fig. 5. — Left and Left Middle: Optical images and the morphological classifications of WIGS galaxies. The image types and filters 
are the same as those listed in Table We present the optical images with two different stretching schemes (inverse hyperbolic sine and 
histogram equalization) that are optimized for demonstrating overall galaxy structure and revealing faint features. Right Middle and Right: 
Velocity maps and velocity dispersion maps derived from the redshifted IFS datacubes. These maps have the same field-of-view in physical 
scale (in parsec) as the optical images. The angular sizes of the velocity and velocity dispersion maps are scaled as if the galaxies are 
observed at 2 : = 1.5. The color bars are in the unit of km s“^. 
























































































11 


High Contrast Image Low Contrast Image Velocity Map Velocity Dispersion Map 




IRAS13120-5453(M4) 
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Fig. 5.— (Continued) 
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IRASF16164-0746(M4) 
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Fig. 5.— (Continued) 
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High Contrast Image Low Contrast Image 
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Fig. 5.— (Continued) 
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IRASF18341-5732 (S) 
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Fig. 5. — (Continued) We note that IRASF18093-5744 “C” refers to the central galaxy located in between the North and South galaxies. 
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Fig. 5.— (Continued) 
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